The low energy spectrum of TeC>2 bolometers: results and perspectives for the 

CUORE-0 and CUORE experiments 
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We collected 19.4 days of data from four 750 g Te02 bolometers, and in three of them we were able 
to set the energy threshold around 3 keV using a new analysis technique. We found a background 
rate ranging from 25 cpd/keV/kg at 3 keV to 2 cpd/keV/kg at 25 keV, and a peak at 4.7 keV. The 
origin of this peak is presently unknown, but its presence is confirmed by a reanalysis of 62.7 kg-days 
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of data from the finished CUORICfNO experiment. Finally, we report the expected sensitivities 
of the CUORE-0 (52 bolometers) and CUORE (988 bolometers) experiments to a WIMP annual 
modulation signal. 

PACS numbers: 07.57.Kp, 29.85.Ca, 95.35.+d 
Keywords: Bolometers, Nuclear Physics, Dark Matter 



I. INTRODUCTION 

Tellurium dioxide bolometers are excellent detectors to 
search for rare processes. Operated at a temperature of 
about 10 mK, these detectors feature an energy resolu- 
tion of a few keV over an energy range extending from a 
few keV up to several MeV. This, together with the low 
level of radioactive background achievable and the low 
cost, makes them ideal detectors for CUORE, an experi- 
ment that will search for neutrinoless double beta decay 
(CVDBD) of 130 Te 0,3. 

CUORE will consist of 988 Te0 2 bolometers of 750 g 
each and is expected to reach a background level at 
the 130 Te Q-value (around 2528 keV [|]|) of the or- 
der of 0.01 counts/keV/kg/y, thus allowing a sensitivity 
to O^DBD down to the inverted hierarchy of neutrino 
masses. CUORE is currently being installed at the Lab- 
oratori Nazionali del Gran Sasso (LNGS) in Italy and is 
scheduled to start taking data in 2014. The experimen- 
tal design was validated by CUORICINO, an array of 62 
Te02 bolometers (for total mass of 40.7 kg) which took 
data at LNGS from 2003 to 2008 and put stringent limits 
on the OidDBD half-life [a, |7J . To test the new assembly 
line and the materials chosen for CUORE, an array of 52 
Te02 bolometers, named CUORE-0, has been recently 
built. It is expected to have a smaller background than 
CUORICINO and will start taking data in Fall 2012. 

Given the high mass, the good energy resolution, and 
the low background, the CUORE-0 and CUORE exper- 
iments can search for rare processes other than OidDBD. 
In particular, provided the energy threshold of few keV, 
they could look also for Dark Matter interactions. Such 
searches were not performed in CUORICINO because the 
energy threshold was too high, of the order of tens of keV. 
A new software trigger was developed to lower the en- 
ergy threshold m. The trigger is based on the matched 
filter algorithm [1, [l(| and also provides a pulse shape 
parameter to suppress false signals generated by detec- 
tor vibrations and electronics noise. 

In this paper we show the energy spectrum of four 
Te02 bolometers originally operated at LNGS in order to 
test the performances of CUORE crystals [ll[ . In three 
of them we were able to set the energy threshold around 
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3 keV, while in the fourth to 10 keV because of the high 
detector noise. The energy spectrum of the three bolome- 
ters exhibits a peak at 4.7 keV whose origin is presently 
unknown. The peak has a constant rate in time, and its 
presence is confirmed by a reanalysis of the data from the 
last two months of operation of CUORICINO. Given the 
observed counting rate, we also evaluate the sensitivity of 
CUORE-0 and CUORE to an annual modulation signal 
induced by WIMP Dark Matter candidates. 



II. EXPERIMENTAL SETUP 

A CUORE bolometer is composed of two main parts: 
a Te02 crystal and a neutron transmutation doped Ger- 
manium (NTD-Ge) thermistor [H E3- The crystal, a 
750 g cube of side length 5 cm, is held by PTFE sup- 
ports in copper frames. The frames are connected to the 
mixing chamber of a dilution refrigerator, which keeps 
the bolometers at a temperature of about 10 mK. The 
thermistor is glued to the crystal and acts as thermome- 
ter. The temperature increase due to an energy depo- 
sition in the crystal is measured by the increase in the 
resistance of the thermistor [l#jl6| . The thermistor is 
biased by a constant current and a change in the volt- 
age across it constitutes the signal [l^j- In the setup 
presented in this paper, called CUORE Crystal Valida- 
tion Run 2 (CCVR2), each crystal was provided with 
two thermistors for redundancy. Two Joule heaters were 
also glued on two crystals and used to inject heat pulses 
with controlled amounts of ener gy t o monitor the detec- 
tor gain and detector efficiencies |18l420l ] . The bolometers 
were operated in the cryogenic R&D facility of CUORE 
whose details can be found in Refs. [2U - I23T ] . 

Data were collected over 23 days with small interrup- 
tions for calibrations and cryostat maintenance. The ef- 
fective live time amounted to 19.4 days. On each of the 
four bolometers, labeled Bl, B2, B3 and B4, we selected 
the thermistor in which the trigger reached the lowest 
energy threshold. 



III. DETECTOR PERFORMANCE 

The energy resolution was evaluated both at the base- 
line level and on the lowest energy peak in the spectrum. 
Randomly triggered events, not containing pulses, were 
used to evaluate the fluctuation of the baseline after the 
application of the matched filter. A Gaussian fit was 
performed on the 30.5 keV Sb line, which is due to EC 
decays of Te isotopes produced by cosmogenic activation 
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during production and air-shipment of the crystals out- 
side the underground laboratory. As one can see from 
Tab. HI baseline and 30.5 keV resolutions are quite close, 
since at these energies the resolution of bolometers is ex- 
pected to be dominated by the noise. The bolometers 
feature excellent energy resolution, with B2, B3, and B4 
below 1 keV FWHM, and Bl at 3.3 keV FWHM. 

The detection efficiencies were measured via an en- 
ergy scan performed at the end of the run by using 
the Joule heaters to provide sequences of pulses from 1 
to 50 keV. However, heater events and actual particle 
events have slightly different pulse shapes which could 
result in slightly different energies and detection efficien- 
cies. The trigger efficiency and any differences that may 
result from the difference between the heater pulses and 
particle events were investigated by Monte Carlo simula- 
tions. 

Heater pulses between 1-50 keV and particle pulses in 
the same energy range were simulated separately using 



the tools described in Ref. 24]. Only B2 and B3 were 



simulated for heater pulses while all the four bolometers 
were considered for the particle pulses simulation. Both 
Monte Carlo runs included simulated detector noise and 
the background from particle pulses, randomly generated 
to match the event rate and the entire measured energy 
spectrum from threshold up to 5.4 MeV, i.e. the energy 
of 210 Po decays in the crystal. We processed the output 
of the simulation in the same way as the heater scan 
measurement. The estimated efficiencies on B2 are shown 
in Fig.[TJ where the good agreement between Monte Carlo 
and data is evident. 

The similarities in the energy dependence demonstrate 
that the trigger acts in the same way on particle and 
heater pulses. This feature is visible also on B3, the 
other bolometer with a heater. We estimate the detec- 
tion efficiencies on the two bolometers without heater (Bl 
and B4) using the Monte Carlo. For each bolometer, we 
evaluate the energy threshold for the data analysis as the 
energy at which the plateau is reached. The detection ef- 
ficiency en was computed as the weighted average of each 
point in the plateau. The energy threshold, the heater- 
measured, and particle-simulated detection efficiencies of 
each bolometer are reported in Tab. |H 



Table I. Baseline resolution (AEb aEC ), resolution at 30.5 keV 
(A.E30 kov), software energy threshold (8e) and detection ef- 
ficiencies measured with the Joule heater ( e ^ eater ) anc j esti- 
mated from the Monte Carlo simulation of the four bolometers 
(4 fC ). 



Bolo 


A_Bbasc | AE30 keV 

[keV FWHM] 


9e 
[keV] 


heater 


MC 


Bl 


3.3 


3.3 ±0.7 


10.0 


n/a 


0.878 ± 0.002 


B2 


0.66 


0.53 ±0.08 


3.0 


0.910 ±0.005 


0.913 ±0.001 


B3 


0.76 


0.83 ±0.09 


2.5 


0.828 ± 0.007 


0.825 ±0.001 


B4 


0.82 


0.59 ±0.12 


2.5 


n/a 


0.828 ± 0.002 
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Figure 1. Detection efficiency on B2: particle Monte Carlo in 
black circles, heater Monte Carlo in gray circles, heater scan 
data in big circles. The shadowed band represents the region 
below the energy threshold chosen for the data analysis. 



IV. ENERGY CALIBRATION 

The energy calibration was performed by inserting two 
wires of thoriated tungsten in proximity of the detector, 
between the cryostat and the lead shields placed exter- 
nally. The main 7 lines of 232 Th, ranging from 511 to 
2615 keV, and the lines generated by metastable Te iso- 
topes in the crystals, ranging from 30 to 150 keV, were 
used for the determination of the calibration function (a 
third order polynomial). The residual difference with re- 
spect to the nominal energy of the peaks in the low energy 
spectrum is shown in Fig. [5J 

The calibration accuracy at energies lower than 30 keV 
was verified by using events from 121 Te and 40 K con- 
taminations in the crystals. These isotopes may decay 
via EC to 121 Sb and 40 Ar, respectively, with the emis- 
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Figure 2. Residuals of the calibration function of B2 on 
metastable Te lines. 
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Figure 3. Summed energy spectrum from B2, B3 and B4 
in time coincidence with an energy deposition on another 
bolometer (including Bl) compatible with the 7 lines from 
40 K (solid) or 121 Te decays (hashed). The dashed lines rep- 
resent the nominal values of the de-excitation energies. No 
event is observed away from these lines. 



Figure 4. Pulse shape indicator (SI) distribution on B2. The 
band at low SI is populated by signal events. At higher values 
there are triggered mechanic vibrations, electronic spikes and 
pileups. The optimal cut for this bolometer was found to be 
SI < 2.2. 



sion of a 7-ray from the daughter nucleus de-excitation 
(507.6 keV or 573.1 keV from 121 Te and 1461 keV from 
40 K) and the de-excitation of the atomic shell Li and K 
{E L1 = 4.6983 keV or E K = 30.4912 keV from 121 Te and 
E K = 3.2060 keV from 40 K). The atomic de-excitation 
process is fully contained in the crystal while, in some 
cases, the 7 can escape and hit another crystal, thus pro- 
ducing a double hit event. We selected double hit events 
requiring that in one of the two crystals there is an energy 
release corresponding to the 7-rays from 121 Te or 40 K, 
and then we looked at the energy of the other crystal in 
the range 0-40 keV (Fig. [3]). Two events from 40 K were 
found (3.04 and 3.18 keV) three from the L x (4.48, 4.67 
and 4.74 keV), and ten from the K de-excitation after 
the 121 Te decay (30.53 ± 0.04 keV). Since all the events 
are compatible with their expected energy, the accuracy 
of the low-energy calibration is confirmed. We also noted 
that the K/L capture ratio for 121 Te is compatible within 
two standard deviations with the expected value of 7 [25[ . 

As a further check of the goodness of our calibration, 
we placed a 55 Fe source next to a bolometer. The X- 
rays produced, with nominal energy between 5.888 and 
6.490 keV, resulted in detected energies that were shifted, 
on average, by only 48 ± 16 eV from their nominal values. 



V. EVENT SELECTION 

We selected signal events using the shape indicator 
variable (SI) described in Ref. 0. This variable is based 
on the x 2 of the fit of the waveforms with the expected 
shape of the signal. The distribution of SI versus energy 
for B2 is shown in Fig.|4] The signal region is easily iden- 
tified, and a considerable part of the pulses generated by 



thermal and electronic noise can be removed with a cut 
on this variable. 

We evaluated the optimal value of the cut on SI us- 
ing the Sb K line at 30.5 keV. We performed a series 
of extended unbinned maximum-likelihood fits, selecting 
the events below a fixed SI value. For each SI cut, the 
spectra of cut accepted events and cut rejected events 
were simultaneously fitted with a Gaussian, representa- 
tive of the signal, plus a first-order polynomial function, 
representative of the background. Signal (es) and back- 
ground (cb) efficiencies after the cut were evaluated di- 
rectly in the fit. For each bolometer we selected the cut 
which maximized the statistical significance, defined as 
£s/\/^Bj corresponding in this case to es = 1. 

VI. MEASURED SPECTRA 

The energy spectra of the four detectors show several 
7 and a peaks that are clearly identified as due to U, Th 
and K contaminations (of the crystals themselves and of 
the experimental setup) and few low energy peaks iden- 
tified as being due to Te metastable isotopes. The region 
between threshold and 40 keV is shown in Fig. [5] In the 
three bolometers with the lowest threshold, a peak ap- 
pears at about 4.7 keV, the energy of the Li atomic shell 
of Sb. This in principle indicates that the line could be 
ascribed to the EC decay of a Te isotope, however none 
of the known or predicted isotope decays can explain our 
observation: 

• The 121 ™Te and its daughter 121 Te decay via EC 
with half-lives of 154 and 17 days, respectively. 
However the intensity of the observed peak is higher 
than the K line of Sb at 30.5 keV. This is in contra- 
diction with the measured value of the K/Li cap- 
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ture ratio for these isotopes which is greater than 
one. 

• The other EC decaying metastable Te isotopes 
(and their daughters) have half-lives smaller than 
4.7 days. As it will demonstrated later, the ob- 
served peak rate is constant over 20 days. 

• 123 Te is a naturally occurring isotope of Tellurium 
(abundance 0.908 ± 0.002 % [26j]) which may decay 
via EC to 123 Sb with a Q-value of 52.2 ± 1.5 keV. 
Since the transition is 2 nd -forbidden unique, it has 
been estimated that very little K capture actually 
occurs and that the majority of EC decays take 
place from the L3 shell. Several searches for Sb K 
lines from EC of 123 Te have been performed. Pos- 
itive evidence, T 1/2 = (1.24 ± 0.10) • 10 13 y, was 
claimed by Watt and Glover [13] but subsequently 
ruled out [H, |2{|. In this paper, we show for the 
first time an energy spectrum down to the L energy 
region. The line we observe however cannot be at- 
tributed to 123 Te, since, as described in the next 
section, the energy is compatible with the Li shell 
(4.6983 keV) and not with the L 3 one (4.1322 keV). 



VII. ANALYSIS OF THE 4.7 KEV PEAK 

In the following we report the analysis of the peak 
in the energy distribution of single bolometers, to pro- 
vide all the possible details and to stimulate a discus- 
sion within the scientific community that will hopefully 
lead to its identification. To determine the intensity, the 
peaking background due to the EC decay of 121m Te- 
121 Te was removed. The 7-rays produced by these iso- 
topes can escape without hitting another crystal, such 
that only the K or the L de-excitations are measured. 
The number of pure L events from 121 Te and 121m Te 
(N121) were estimated using a Monte Carlo simulation 
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Figure 6. Best fits at the 4.7 keV line on the three bolometers 
with low threshold. On the bottom right, the combined profile 
negative log-likelihood projected over the number of signal 
events. 



Table II. Best fit results for the 4.7 keV line and estimated 
peaking background from 121 Te. The error on the energy 
includes the systematic due to the calibration. 



Bolometer 


Energy 


FWHM 


Nsig 


Nm 




[keV] 


[keV] 


[counts] 


[counts] 


B2 


4.75 ± 0.28 


0.57 ±0.11 


191±| 


6.1 ±1.3 


B3 


4.66 ±0.28 


0.87 ±0.10 




13.5 ± 1.6 


B4 


4.76 ± 0.38 


0.75 ±0.12 


205t« 


7.5 ±1.4 



Figure 5. Low energy spectra of the four bolometers. Full 
statistics (19.4 days) with pulse shape cut applied. 



based on GEANT4 30 . normalizing the single hit spec- 
tra to the measured rate of the most intense 121m Te peak 
(294.0 keV). 

To estimate the intensity and the energy of the line 
from the data, we performed a separate extended un- 
binned maximum-likelihood fit for each bolometer, using 
a likelihood function constituted by a Gaussian plus two 
exponential functions to reproduce the background. We 
set the pulse shape cuts at the estimated optimal val- 
ues. The selection efficiencies were recomputed on the 
peak at 4.7 keV, since the events in this peak are more 
plentiful than in the Sb K peak at 30.5 keV. The ob- 
tained number of events N e was then corrected taking 
into account the detection and cut efficiencies and the 
expected background A421, according to the equation 
N sig = (N e — N121)/ (eo^s)- Best fits are shown in Fig. [6] 
In Tab. [II] we report the summary of the peak parameters 
of each bolometer, also including the estimated peaking 
background N\2i- 

We combined the profile negative log-likelihoods of N e 
and iVi2i of the three bolometers (Fig. [5]) to compute 
the average rate. The estimated value of the number of 
signal events is 223 ± 22 counts/crystal, from which we 
evaluated the line intensity in Te02 to be I = 15.3 ± 
1.5 counts/day /kg. The energy averaged over the three 
bolometers is 4.72 ± 0.18 keV. 
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Figure [7] shows the peak rate along the duration of the 
data-taking, which is in a very good agreement with a 
constant distribution, indicating that this signal is not 
due to a short-living radioactive contamination. It also 
indicates that any variation of the detection efficiency 
with time is negligible compared to the statistical error. 
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Figure 8. Comparison of CCVR2 (black circles) and CUORI- 
CINO (red triangles) data, with a zoom on the 4.7 keV. Un- 
like in CCVR2, CUORICINO data are not corrected for the 
efficiencies and the calibration at low energies is not opti- 
mized. 



Figure 7. Rate of the 4.7 keV peak versus solar time in all 
three CCVR2 bolometers summed together. No variation is 
appreciable. 

To check the stability of the line over a longer period of 
time, we checked its presence in the CUORICINO data. 
During the last two months of operation of CUORICINO, 
the data-acquisition system being developed for CUORE 
(the same used for CCVR2), was run in parallel with 
the old system, saving to disk the continuous stream 
of data. The live-time of continuous data amounts to 
47.5 days. We ran the new trigger on these data but we 
were able to reach a threshold below 4 keV only on 4 
bolometers, because of the high vibrational noise trans- 
mitted by the holder to the crystals (CCVR, CUORE-0 
and CUORE holders have been explicitly redesigned to 
lower the transmitted noise). These 4 bolometers had 
smaller size (3x3x6 cm 3 , 330 g) than other CUORICINO 
and CUORE bolometers, and featured higher signal to 
noise ratio. 

The sum energy spectrum of the four bolometers is 
overlaid to the sum energy spectrum of the CCVR2 
bolometers in Fig. [5] In CUORICINO the heater energy 
scan was not performed, therefore the spectrum shown is 
not corrected for the detection efficiencies. As it can be 
seen from the figure, the peak is very well visible, with a 
fitted energy resolution of 0.52 ± 0.04 keV FWHM. The 
intensity and the average energy are found to be similar 
to the CCVR2 ones, being (10.0 ± 0.6)/e counts/day/kg 
and 4.73 ± 0.01 keV, respectively. We also note that the 
30.5 keV peak disappeared, because the 121m - 121 Xe iso- 
topes decayed away during the 5 years of underground 
data taking. For this reason, the calibration is not as 
accurate at low energies. However, the error is expected 
to be less than 0.5 keV. 



VIII. SENSITIVITY TO WIMPS 

It has to be remarked that Te02 bolometers, unlike 
other bolometric detectors for Dark Matter, do not dis- 
criminate the nuclear recoils induced by WIMPs from 
the radioactive /3/7 background. Nevertheless, the high 
mass and low background achievable with these detectors 
make it possible to search for an annual modulation of 
the counting rate. 

Compared to CUORICINO, the CCVR2 rate has the 
same behavior at energies greater than 10 keV, but is 
considerably higher at lower energies (Fig. [5J. We are 
unable to explain this difference at present, however we 
expect that the CUORE-0 and CUORE low energy back- 
ground will not be higher than the CCVR2 one. All the 
materials used in detector construction, in fact, will be 
the same as those employed for CCVR2. Moreover, in 
the case of CUORE-0, the bolometers will be operated 
in the same cryostat of CUORICINO, which has lower 
radioactive contaminations compared to the one used to 
operate CCVR2. To be conservative, we estimate the 
CUORE-0 and CUORE sensitivity to WIMPs assuming 
the background rate of these experiments to be equal to 
the one measured on CCVR2. We assume that the noise 
will be under control and that all bolometers will achieve 
a 3 keV threshold. We focus on the energy region be- 
tween threshold and 25 keV, featuring an observed back- 
ground counting rate ranging from about 25 cpd/keV/kg 
to 2 cpd/keV/kg. 

We perform toy Monte Carlo simulations generating 
background events from the CCVR2 fit function shown 
in Fig. and WIMP events from the predicted distri- 
bution described in Ref. [3l|, using the following WIMP 
parameters: density pw — 0.3 GeV/cm 3 , average veloc- 
ity vq — 220 km/s and escape velocity from the Galaxy 
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Figure 9. WIMP signal: difference between the 3-month in- 
tegrated spectra across December the 2nd and June the 2nd, 
for a WIMP cross section of 1CP 41 cm 2 and masses between 
5 and 100 GeV. 



v esc — 600 km/s. The quenching factor (QF) of the in- 
teractions in TeC>2 is set to 1 [32j . We include the depen- 
dence of the WIMP interaction rate on the time in the 
year, and estimated the background+signal asymmetry 
subtracting the 3-month integrated spectrum across De- 
cember the 2nd from the 3-month integrated spectrum 
across June the 2nd. The resulting differential spectrum 
is fitted with the expected shape induced by the mod- 
ulation, H\ (examples are given in Fig. [H]), and with a 
flat line at zero counts, Hq. The cross section in the toy 
simulation is lowered as long as, in a set of experiments, 
the fit probability of the Hi hypothesis is greater than 
the H one at least 90% of the times. 

This procedure defines the cross section that could be 
sensed for a fixed WIMP mass. The 90% CL upper limit 
sensitivity to the cross-section as a function of the WIMP 
mass is reported in Fig. [10] for 3-years of CUORE-0 data- 
taking and 5-years of CUORE. The comparison with 
other experiments shows that CUORE-0 could test the 
indication of a - 10 GeV WIMP from the DAM A (no- 
channeling), CoGeNT and CRESST experiments, while 
CUORE could completely test the DAMA results, under 
the hypothesis that Dark Matter is made of WIMPs. We 
reiterate that because the quenching factor for nuclear 
recoils compared to electron recoils in Te02 bolometers 
is 1, the 2-6 keV energy region of DAMA corresponds 
to 7-20 keV assuming scattering on Na (QF=0.3) or 22- 
67 keV assuming scattering on I (QF=0.9). Therefore, 
in Te02 bolometers it will be possible to look at lower 
energies and to study with larger detail the shape of the 
modulation spectrum, thus providing new information to 
this complicated search. 




WIMP mass [GeV/c 2 ] 

Figure 10. 90% sensitivity to WIMPs of CUORE-0 and 
CUORE assuming a 3 keV threshold for all detectors and the 
same background level of the CCVR2 detectors. Evidences of 
DAMA-3cr [13, CoGeNT-90% and CRESST-2cr [H are 
reported for comparison. 
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